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Based on high-order harmonic generation (HHG) spectra obtained from solving the time- 
dependent Schrodinger equation for atoms, we established quantitatively that the HHG yield can be 
expressed as the product of a returning electron wave packet and the photo-recombination cross sec- 
tions, and the shape of the returning wave packet is shown to be largely independent of the species. 
By comparing the HHG spectra generated from different targets under identical laser pulses, accu- 
rate structural information, including the phase of the recombination amplitude, can be retrieved. 
This result opens up the possibility of studying the target structure of complex systems, including 
their time evolution, from the HHG spectra generated by short laser pulses. 
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When an atom is subjected to a strong driving laser 
field, one of the most important nonlinear response pro- 
cesses is the generation of high-order harmonics. In the 
past decade, high-order harmonic generation (HHG) has 
been used for the production of single attosecond pulses 
P, 0, 0] and attosecond pulse trains thus opening 
up new opportunities for attosecond time-resolved spec- 
troscopy. HHG is understood using the three-step model 
(TSM) d @, 01 - first the electron is released by tun- 
nel ionization; second, it is accelerated by the oscillating 
electric field of the laser and later driven back to the tar- 
get ion; and third, the electron recombines with the ion 
to emit a high energy photon. A semiclassical formula- 
tion of the TSM based on the strong-field approximation 
(SFA) is given by Lewenstein et aZ [7[. In this model (of- 
ten called Lewenstein model), the liberated continuum 
electron experiences the full effect from the laser field, 
but not from the ion that it has left behind. In spite 
of this limitation, the SFA model has been used quite 
successfully, in particular, for analysis of the attosecond 
synchronization of high harmonics, see Mairesse et al Q 
and references therein. However, since the continuum 
electron recombines when it is near the parent ion, the 
neglect of electron-ion interaction in the SFA model is 
rather questionable. 

According to the TSM, the last step of HHG is analo- 
gous to the radiative recombination process in electron- 
ion collisions. Thus one may write the HHG signal as 



S{uj) = W{E) X \d(uj)\^ 



(1) 



where d{Lo) is the photo-recombination (PR) transition 
dipole and W{E) is the returning "electron wave packet" . 
Electron energy E is related to the emitted photon energy 
oj hy E = oj — Ip, with Ip being the ionization potential 
of the target. Clearly the HHG signal S{uj) and W{E) 
depend on the laser properties. On the other hand, d{uj) 



is the property of the target only. The factorization in 
Eq. (1) is most useful when one compares the HHG spec- 
tra from two different targets in the identical laser field. 
Assuming that the shape of W{E) is species indepen- 
dent, by measuring the relative HHG yields, one can de- 
duce the PR cross section of one species if the PR cross 
section of the other is known. The validity of Eq. (1) 
has been shown recently in Morishita et al. |9i] using 
HHG spectra calculated by solving the time-dependent 
Schrodinger equation (TDSE) for atoms. The validity of 
this factorization has also been shown for rare gas atoms 
by Levesque et al. and for N2 and O2 molecules 11 1 , 
where the HHG spectra were calculated using the SFA 
model. In the SFA the continuum electron is approxi- 
mated by plane waves, thus the dipole matrix elements 
are calculated in the plane wave approximation (PWA). 

In this Letter, we have two goals. The first is to show 
that electron wave packets obtained from the SFA model 
and from the TDSE calculation are nearly identical, but 
the transition dipoles calculated from PWA are signifi- 
cantly different from using scattering waves (SW). This 
result suggests a scattering wave based strong-field ap- 
proximation (SW-SFA) for harmonic generation where 
the wave packet is derived from the SFA but the transi- 
tion dipole is calculated using accurate SW. The second 
goal is to check whether Eq. (1) can be extended to the 
level of complex amplitudes such that one can relate the 
phases of HHG to the phases in the transition dip oles. 
Since phases of harmonics can be measured ij, [l3| , and 
they are needed in order to incorporate the effect of prop- 
agation in the macroscopic medium, such a study is im- 
portant. 

First in Fig. 1 we compare the PR cross sections of 
Ar, Xe and Ne calculated by treating the continuum 
electrons using PWA to results calculated with accurate 
SW's. Clearly they show significant differences. They 
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FIG. 1: (Color online) Photo-recombination cross sections of 
Ar (a), Xe (b) and Ne (c), obtained by using exact scatter- 
ing wavefunctions (solid curves) and within the plane-wave 
approximation (dashed curves) for the continuum electrons. 
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FIG. 2: (Color online) Comparison of the returning electron 
"wave packets" extracted from the HHG spectra of Ne, ob- 
tained by solving the TDSE (solid black line), and from the 
SFA model (dashed blue line). Also shown is the TDSE result 
for the wave packet from scaled H (dotted red line). For laser 
parameters, see text. 



reflect the well-known facts that plane waves are poor 
approximations for representing continuum electrons in 
atoms and molecules for energies in the energy range of 
tens to hundreds eV's. 

To check whether the target structure affects the re- 
turning electron wave packet, in Fig. 2 we compare the 
wave packets W{E) for Ne deduced from the TDSE and 
SFA results, using Eq. (1). In the SFA case, the transi- 
tion dipole is calculated within the PWA. Also shown is 
the W{E) obtained from scaled atomic hydrogen, with 
the effective nuclear charge chosen such that the ioniza- 
tion potential of its Is ground state is the same as of 
Ne(2p). We used a laser pulse with duration (FWHM) 
of 10.3 fs, peak intensity of 2 x 10^"* W/cm^ and mean 
wavelength of 1064 nm. Note that we have normalized 
the results near the cut-off. The normalization is to ac- 
count for the difference in the tunneling ionization rates 
from SFA and from TDSE, or from the different species. 
This comparison shows that the shape, or the energy de- 
pendence of the electron wave packets, depends only on 
the laser parameters. 

Having established that the wave packet can be ob- 
tained from the SFA model, we now examine the accu- 
racy of HHG calculated using the SW-SFA model where 
the wave packet is extracted from the SFA model and the 
transition dipoles are calculated using SW. In Fig. 3 we 
show the HHG spectra obtained from the TDSE, SFA, 
and SW-SFA for Ar, Xe, and Ne. For Ar and Ne, the laser 
pulse has peak intensity of 2 x 10^* W/cm^ and mean 
wavelength of 800 nm. The laser duration (FWHM) is 
10 fs for Ar and 20 fs for Ne. For Xe, the corresponding 
parameters are 5 x 10^^ W/cm^, 1600 nm, and 7.8 fs, re- 
spectively. The HHG yields for Ar are shifted vertically 
in order to show their detailed structures. For Ne and 
Xe, the SFA and SW-SFA results are normalized to the 
TDSE results near the cutoff, i.e., close to 3.2Up + Ip, 



where Up is the ponderomotive energy. 

The results in Fig. 3 clearly demonstrate the good 
improvement of the SW-SFA over the SFA in achiev- 
ing better agreement with the TDSE results. Since the 
SFA gives the correct wave packet, its prediction would 
be "reasonable" in the energy region where the dipole 
matrix element is rather flat, i.e., in the higher photon 
energy region. Thus the SFA would give adequate pre- 
diction of the HHG spectra usually near the cutoff re- 
gion (after spectra are renormalized). This fact has been 
known ^. The improvement of SW-SFA occurs usually 
at lower photon energies where the PWA for the con- 
tinuum electron is grossly incorrect. In particular, the 
transition dipole from PWA goes through zero at some 
lower energies, see Fig. 1. This is the energy region where 
the SFA suffers the largest errors. Because of the zeros 
in the dipole matrix elements in the PWA, the deduced 
wave packets from SFA would suffer large errors at the 
corresponding energies. These errors are reflected as the 
sharp spikes in the HHG spectra calculated using the 
SW-SFA model. 

For a realistic description of the experimental harmonic 
spectra, the effect of phase matching and macroscopic 
propagation should be addressed. To this end, the knowl- 
edge of the harmonic phase is necessary. Thus extending 
Eq. (1) to include the phase, can the phase of the har- 
monic be expressed as the sum of the phase from the wave 
packet and from the PR transition dipole? First wc es- 
tablish that there is a close relationship between the har- 
monics phase (f> and the PR dipole phase 5. To be specific, 
we focus on Ar target. We calculated the phase differ- 
ence A(j) for each harmonic generated from Ar and from 
its scaled hydrogen (reference) partner under the same 
laser pulse. These calculations were carried out using 
the TDSE with 4-cycle and 10-cycle laser pulses, intensi- 
ties of 1 and 2 x lO^"' W/cm^, and wavelength of 1064 nm 
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FIG. 3: (Color online) Comparison of the HHG yields ob- 
tained from numerical solution of the TDSE (solid red lines), 
the SFA (dashed blue lines), and the SW-SFA model (solid 
black lines) for Ar (a), Xe (b), and Ne (c). For laser parame- 
ters, see text. 



and 800 nm. In Fig. 4(a) we compare Acj) — (f>^^ ~ (f>^'^-f 
with the PR dipole phase difference AS = S^^ — 5^'^^ . 
Here we have shifted the harmonic phase difference to 
match the PR dipole phase difference at = 60 eV. 
Clearly, the two agree very well for the different lasers 
used. In particular, the phase jump near 40 eV (due to 
the Cooper minimum in Ar) is well reproduced. This 
indicates that the phase of the wave packets from the 
two systems are almost identical (up to a constant shift). 
Similar agreements were also found for Xe and Ne, as 
shown in Fig. 4(b) and (c), respectively. Here the laser 
pulses of 4 cycles duration are used, other parameters are 
given as shown in the labels. This result allows one to 
obtain the harmonic phase (j) from the harmonic phase of 
the partner atom cff'^^ by using (f) — (ff^^ + A5. 

We have also applied the same procedure by compar- 
ing the TDSE and SFA results for the same target and 
found that A0 — cj)^^^^ — cf)^^^ no longer agrees well 
with AS = - S^^ . This indicates that the phase of 
the electron wave packet calculated from SFA differs from 
the one calculated by TDSE, even though their magni- 
tudes agree well. How significant these differences affects 
the HHG spectra after macroscopic propagation? To this 
end we calculate the HHG spectra by coherently averag- 
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FIG. 4: (Color online) Extracted harmonic phase difference 
Acj) between Ar and scaled hydrogen obtained with different 
lasers as function of emitted photon energy. The PR dipole 
phase difference A(5 is given as solid black line, (a) Ar, (b) 
Xe,and (c)Ne. Io=10" W/cm^ 



ing the induced polarization over an intensity range of 
the driving laser. In Fig. (5) we show the results for Ar 
from the TDSE, the SW-SFA, and the one with the wave 
packet extracted from the scaled hydrogen. All of these 
results are coherently averaged over 11 equally-spaced in- 
tensities in the range from 1.8 to 2.2 x IQ-^'' W/cm^. The 
laser is of 800 nm wavelength and 30 fs (FWHM). The 
scaled H result is indeed in quite good agreement with 
the exact TDSE calculations. This is not surprising since 
we have shown that the phase of the wave packet from 
the scaled H and from Ar are almost identical at a single 
intensity. For SW-SFA, the agreement is not as good, 
but the improvement over SFA is still significant. The 
phase in SFA (or SW-SFA) can probably be improved by 
adding some correction to the semi-classical action, for 
example, as has been suggested [13; At present, it is 
better to extract the phase of the wave packet from the 
companion atomic target where TDSE calculations can 
be carried out. 

Here we comment on the computational details. The 
solution of the TDSE and the choice of one-electron 
model potential for describing the atom have been de- 
scribed previously 3, 



1, 



[ . The electric field of the laser 
pulse is written in the form E{t) = EQa{t) cos{u!t) , with 
the envelope given by a{t) = cos^(7rt/r), where t is 2.75 
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FIG. 5: (Color online) HHG spectra for Ar from the "sim- 
ulated" macroscopic propagation. Shown are results from 
the exact TDSE (solid red line), SW-SFA(dotted black), and 
by using the wave packet extracted from TDSE solution for 
scaled H(ls) (dashed blue). For laser parameters, see text. 



times the FWHM of the laser pulse. To calculate the PR 
cross section, the scattering wavefunction is expanded in 
terms of partial waves and the transition dipole is calcu- 
lated for the continuum electron that has the wave vector 
along the polarization axis only. 

Before concluding we mention several earlier related 
works. There exists a wealth of literature aiming at im- 
proving the SEA model, e.g., by including Coulomb dis- 
or by eikonal approximations [iF 
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18j . or by eikonal approximations [15[. In 
these approaches, the PR processes are still treated ap- 
proximately. For example, use of Coulomb wave for the 
continuum electron would not produce the Cooper min- 
imum in the PR cross section in Ar (see Fig. 1). The 
advantage of SW-SFA is that it factors out the target 
structure explicitly. A minimum in the HHG spectra may 
be attributed to the minimum in the PR cross section 
and this position should not change with laser parame- 
ters. Such minima are of particular interest for molecu- 
lar targets since minima in the molecular dipole matrix 
element may be interpreted as due to the interference 
between the emission amplitudes from different atomic 
centers. Interference minima have been observed exper- 
imentally in CO2 by different groups [l^, [l^], but the 
observed positions of the minimum are not identical and 
thus other possible interpretations have been suggested 
2l| . Another hot topics in recent years is the tomo- 
graphic method for imaging the molecular orbitals [23 |. 
This pioneering work deduced the dipole matrix elements 
of N2 molecules by comparing the HHG spectra of N2 
vs Ar, using the factorization Eq. (1). In order to use 
the tomographic procedure to obtain the ground state 
wavefunction of N2, they approximated the continuum 
wavefunctions of Ar and N2 by plane waves. In view of 
Fig. 1, their success of extracting good valence orbital 
wavefunction of N2 is surprising. We note, however, in 
Itatani et al [22| . the continuum electron energy is set 
equal to the photon energy, arguing that the electron re- 



combining near the core should gain the additional bind- 
ing energy. For Ar, this would shift the PWA curve in 
Fig. 1 by 15.7 eV, making the PWA result much closer to 
the SW result. However, this shift does not always work, 
see Xe and Ne examples in Fig. 1. 

In conclusion, we have established quantitatively that 
the last step of the three-step model of HHG can indeed 
be expressed as the photo-recombination (PR) process of 
the returning electron wave packet. The wave packet de- 
pends nonlinearly on the laser, but its shape and phase 
are largely independent of the target. Thus if the PR of a 
reference target is known, the PR of another target can be 
derived by measuring the HHG of the two species under 
identical laser pulses. Since the results should be inde- 
pendent of the lasers, this allows for an important check 
on the accuracy of the measurements. We also showed 
that the HHG spectra can be calculated using the SW- 
SFA model. This model describes well the single atom 
HHG intensity, but the phase needs further corrections. 
For complex systems, SW-SFA would be a good start- 
ing point for describing the HHG spectra since the PR 
process is accurately incorporated. While our conclusion 
has been derived based on atomic targets and in the sin- 
gle active electron model, we anticipate that the results 
are applicable to molecules where accurate TDSE calcu- 
lations are not available in general. The present result of- 
fers a systematic roadmap for extracting target structure 
information from the high-order harmonics generated by 
intense lasers. 
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